High-order harmonic generation with resonant core excitation by ultraintense x rays 
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High-order harmonic generation (HHG) is combined with resonant x-ray excitation of a core 
electron into the transient valence vacancy that is created in the course of the HHG process. To 
describe this setting, I develop a two-active-electron quantum theory for a single atom assuming no 
Coulomb interaction among the electrons; one electron performs a typical HHG three-step process 
whereas, another electron is excited (or even Rabi flops) due to the intense x rays from the core 
into the valence hole and after the first electron has left the atom. Depending on the probability to 
find a vacancy in the valence and the core, the returning continuum electron recombines with the 
valence and the core emitting HHG radiation that is characteristic for the combined process. After 
presenting the theory of x-ray boosted HHG for continuous-wave light fields, I develop a description 
for x-ray pulses with a time varying amplitude and phase. My prediction offers novel prospects 
for nonlinear x-ray physics, attosecond x rays, and HHG-based time-dependent chemical imaging 
involving core orbitals. 

PACS numbers: 42.65.Ky, 32.80.Aa, 32.30.Rj, 41.60.Cr 



I. INTRODUCTION 

Since the surprise discovery of high-order harmonic 
generation (HHG) by atoms in intense optical laser 
fields p], [|| , HHG has spawned the field of attoscience, 
is used for spectroscopy, and serves as a light source in 
many optical laboratories 0-0|- Initially, there has been 
a period of learning about the fundamental physics of 
HHG which lead to a phenomenological description of 
HHG by a single atom in terms of a three-step model 
and quantum theories • Soon it was realized that 

the high harmonic (HH) light from a dense gas is not 
well-described by the HH spectrum of a single atom and, 
additionally, one needs to account for propagation ef- 
fects of the optical laser and the HH radiation in the 
macroscopic medium pil-23]. Nowadays, more and more 
applications based on HHG have come into focus [sec, 
e.g., Refs. H, H, and0j- To name a few, in atoms, there 
is the phase measurement of resonant two-photon ion- 
ization in helium (24j and intense x-ray generation [23j j ; 
for molecules, there is tomographic imaging of orbitals 
and related methods for structure determination [25H30j , 
and the control of electron localization in photodissoci- 
ation 31]. Finally, for condensed matter, an attosecond 
spectroscopy has been developed 133. 

Present-day theory of HHG [1, H 0] mostly gravitates 
around HHG from valence electrons and the single-active 
electron (SAE) approximation (HHl]. Within the SAE 
approximation and the three-step model of HHG [1, [t| , 
the optical laser tunnel ionizes a valence electron and ac- 
celerates it in the continuum. When the optical laser 
field changes direction, the liberated electron is driven 
back to rescatter with the parent ion. This may cause 
the electron to recombine with the ion whereby the ex- 
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cess energy due to the atomic potential and due to the 
energy gained from the optical laser field is released in 
terms of a high harmonic of the optical frequency. In 
HHG each rescattering generates attosecond bursts of ra- 
diation. HHG from an ultrashort optical laser pulse leads 
to an attosecond pulse train [35| in which the individual 
attosecond pulses are separated by a half cycle period 
of the optical laser from which a single attosecond burst 
can be isolated by filtering out the harmonics close to the 
HHG cutoff (maximum photon energy for a given optical 
laser intensity) [36| . By optimizing the quantum path of 
the continuum electron, the HHG cutoff can be increased 
by a factor of 2.5 $t\ . 

HHG beyond the SAE model with only a single atomic 
valence state has been studied in a few works. In 
Refs. [38H4ll . HHG from a coherent superposition of the 
ground state and an excited state of an atom was investi- 
gated. A shift of the HH spectrum was found amounting 
to the energy difference between the two states involved 
with the shifted plateau having a similar height as the one 
they get for HHG from an electron either in the higher 
state or the lower energy state alone. Recently, it was 
both experimentally observed and theoretically explained 
in [42|, |43j that molecular states below the highest occu- 
pied molecular orbitals have significant impact on the 
HHG emission; the authors of [44j investigate in strong- 
field approximation pH l4ol |46| the influence of quantum 
interference on HHG spectra of molecules with two opti- 
cally active electrons from different states. Multiple cut- 
offs from plasmon-like excitations in ions were predicted 
in Ref. 47t The role of multielectron effects regarding the 
recombination probability of the recolliding electron was 
discussed in Ref. H3 and HU. A two-electron scheme was 
considered that uses sequential double ionization by an 
optical laser with a subsequent nonsequential double re- 
combination; in helium it leads to a second plateau with 
about 12 orders of magnitude lower yield than the pri- 
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mary HHG plateau [49| . 

Frequently, the SAE approximation and valence- 
electron HHG is also applied to two-color HHG where 
a optical laser is combined with vuv/xuv light from an- 
other HHG process [5(1 l5l| : the vuv radiation assists, 
thereby, in the ionization process leading to an overall 
increased yield [H(l IH, [HI . This principle is evolved fur- 
ther by using attosecond xuv pulses to assist in the HHG 
process which increases the HHG yield for a certain fre- 
quency range by enhancing the contribution from specific 
quantum orbits [54|-|59l|. H owever, there are only few ex- 
ceptions, e.g., Refs. l27l. El. l60l - [63l in which many-electron 
effects are treated for two-color HHG. In Refs. [(3(3 andlr3ll 
the influence of other electrons is included implicitly by 
using a frequency-dependent polarizability for the atoms; 
the xuv light is found to cause new plateaus to emerge at 
higher energies, however, with a much lower HHG yield. 

Alternatively to HHG sources for xuv light and, par- 
ticularly, x rays, there are synchrotrons and — most rel- 
evant for this work — the newly constructed free elec- 
tron lasers (FEL) that provide exciting novel opportu- 
nities for strong- field physics |64| . For example, the 
Free Electron Laser in Hamburg (FLASH) [65j offers 
xuv to soft x-ray radiation that is ideal for studying 
the valence and shallow core of atoms or the Linac 
Coherent Light Source (LCLS) [H, [(37j which delivers 
soft and hard x rays suitable for x-ray diffraction with 
atomic resolution and examining core electrons of most 
elements in the periodic table. Existing FELs produce 
light of unprecedented intensity which can be used — 
among many other applications — to manipulate optical 
strong-field processes. They operate according to the 
self-amplification of spontaneous emission (SASE) FEL 
which produce chaotic light [68l - [70| that can be modeled 
by the partial coherence method (PCM) (Zlrlzl- 

The combination of optica l light and x rays has proven 
very beneficial in the past [7^14 7a |. e.g., to control the 
interaction of the x rays with optical light in various 
ways involving atoms and molecules (please see Refs.l77l- 
l83l and Refs. therein). Recently, x-ray and optical 
wave mixing — specifically sum-frequency generation of 
x rays and optical light — was demonstrated at LCLS by 
Glover et al. [HI with a conversion efficiency of 3 x 1CP 7 
which offers a novel way to probe matter. Such x-ray 
and optical wave mixing was proposed theoretically [85h - 
|87| many years ago in the 1970s but only recently it 
has been observed experimentally due to the construc- 
tion of FELs with unprecedented x-ray intensity [84j . 
The experiment of Glover et al. is a very encourag- 
ing motivation of theoretical research on x-ray boosted 
HHG [il, [H HI] . Within the SAE, HHG from core elec- 
trons of neon that are ionized by x rays was studied the- 
oretically [S8j | . A single- atom efficiency comparable with 
conventional HHG was found and attosecond pulses were 
isolated from the x-ray-boosted HHG spectra. 

In this work, I would like to develop in detail the 
theory of HHG by an optical laser with an additional 
resonant excitation (or even Rabi flopping) in intense 



x rays. Specifically, I will focus on the case of continuous- 
wave (CW) light, first, and a CW optical laser in com- 
bination with arbitrary x-ray pulses, second. The prin- 
cipal idea is sketched in Fig. [1] In the parlance of the 
three-step model [1, Q , the x-ray-boosted HHG process 
proceeds as follows: first, the atomic valence is tunnel 
ionized; second, the liberated electron propagates freely 
in the electric field of the optical laser; third, the direc- 
tion of the optical laser field is reversed and the elec- 
tron is driven back to the ion and eventually recombines 
with it emitting HHG radiation. The excursion time of 
the electron from the ion is approximately 1 fs for typ- 
ical 800 ran laser light. During this time, one can ma- 
nipulate the ion such that the returning electron sees 
the altered ion as depicted in Fig. Q] Then, the emit- 
ted HHG radiation bears the signature of the change. 
Perfectly suited for this modification during the propa- 
gation step is x-ray excitation of an inner-shell electron 
to the valence shell. The recombination of the return- 
ing electron with the core hole leads to a large increase 
of the energy of the emitted HHG light as the energy of 
the x-ray photons is added to the unmodified HHG spec- 
trum. A prerequisite for this to work certainly is that 
the core hole is not too short lived [64|, i.e., it should 
not decay before the electron returns. In Ref. [62, we 
examined HHG in krypton atoms where tunnel ioniza- 
tion leads to 4p vacancies. The xuv light is tuned to 
the 3d —> 4p resonance in the cation which leads to a 
second high-yield HHG plateau that is shifted to higher 
energies by the xuv photon energy. Similarly, we studied 
x-ray-boosted HHG for neon 2p vacancies where x ray s 
were tuned to the Is — > 2p resonance in the cati on J63B . 
If the x rays are very intense, even Rabi flopping |89h93| | 
of the core electron is possible. Such Rabi oscillations 
have been investigated in neon atoms (without optical 
laser) (^[1111]. 

Our scheme offers new prospects for HHG involving 
core electrons. However, if one is only interested in an 
extension of the HHG cutoff into the kiloelectronvolt 
regime, one may use conventional HHG that is labora- 
tory size and valence-electron-based. In this case, the 
optical laser intensity needs to be high ~10 15 W/cm 2 and 
its wavelength should be in the midinfrared instead of the 
nir. The longer wavelength, however, reduces HHG ef- 
ficiency due to longer continuum propagation and thus 
enlarged spreading of the wave packet. Additionally, 
such a high optical las er in tensity increases the ionization 
rate substantially [§2, Il00j | leading to large electron back- 
ground that causes phase mismatching. A good HHG 
yield can, nonetheless, be obtain ed by a judicious choice 

of the gas pressure [23. IToil4To3 | . 

The paper is structured as follows. The theory of HHG 
in the presence of a CW optical laser and CE x rays is 
developed in Sec. [TT| I discuss the wavefunction ansatz 
[Sec. Ill A| . the Hamiltonian [Sec. Ill B| . the equations of 
motion (EOMs) [Sec. Ill C| . the time-dependent dipole 
transition matrix element [Sec. HID] , and the HHG spec- 
trum [Sec. Ill E] , A formalism for a CW optical laser and 
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(a) 



(b) 




FIG. 1. (Color online) Schematic of the three-step model of 
HHG [1, [i| augmented by x-ray excitation of core electrons: 
(a) tunnel ionization of a valence electron in the optical laser 
field; (b) while the first electron propagates freely in the con- 
tinuum, a second electron from a core orbital is excited to 
the valence hole; (c) eventually, the first electron is driven 
back to the residual ion by a reversal of the direction of the 
optical laser field and recombines with the newly formed core 
vacancy emitting HHG light. 



arbitrary x-ray pulses is derived in Sec. IIIII I present the 
Hamiltonian, the wave function ansatz, and the EOMs 
[Sec. IIII A] , decouple and integrate the EOMs [Sec. lfflB"] . 
and calculated the time-dependent electric dipole mo- 
ment and the HHG spectrum [Sec. IIII Cj . Finally, con- 
clusions are drawn in Sec. IIV1 In the Appendix, I dis- 
cuss the harmonic photon number spectru m (H PNS). All 
equations are formulated in atomic units [1041 ] . 



II. HIGH HARMONIC GENERATION WITH 
CONTINUOUS WAVE X RAYS 

I consider the HHG process for an atom in the two- 
color light of an optical laser and intense x rays. The 
optical laser-only problem was treated by Lewenstein et 
al. [ll| for a valence electron in SAE [H H]. Here 
I consider also a core electron in addition to the va- 
lence electron [Fig.Q]. The combined treatment of HHG 
and x rays, which induce Rabi flopping [73|, I89U98I ] in 
the residual ion, represents a genuine two-active-electron 
problem. 



A. Symmetry and basis states 

The theory is developed here for a closed-shell atom 
with a spin-singlet ground state. I set out from a non- 
rclativistic, independent-electron solution of t he atomic 
electronic structure of Hartree-Fock-Slater type [l05lll06j | 
with spherical averaging of the one-electron potential 
leading to a cent ral p otential as described by Herman 
and Skillman (77l . Il07j |. Thus, for an isolated atom with- 
out light fields, atomic orbitals with the same principle 
and angular quantum numbers are degen erate with re- 
spect to the magnetic quantum number [l08j |. As no 
spin-dependent terms are considered in the Hamiltonian 
for the atom in the optical and x-ray light, all quanti- 
ties can be expressed in terms of spatial atomic orbitals. 
Observables, however, need to be multiplied by a factor 
of 2 in order to account for the two electrons per spatial 



orbital. 

The essential states relevant to describe HHG in two- 
color light are deduced as follows. I assume that the po- 
larization vectors of the optical laser e\, and the x rays ex 
are linearly polarized along the z axis e z , i.e., e\ = ex = 
e z holds. In this case, the magnetic quantum number is 
conserved in the interactions in the absence of spin-orbit 
coupling whic h influenc es strong-field ionization of heav- 
ier atoms [78l Il09l - llll |. However, strong optical fields 
align the atomic orbitals along the linear optical laser 
polarization axis. I will not treat spin-orbit coupling 
in what follows and assume that the magnetic quan- 
tum number is a good quantum number. Hence, one- 
electron states of relevance to the problem are the valence 
states | a; m) with magnetic quantum numbers to e M a 
with M. a = {— l a , . . . , l a } for principal n a quantum num- 
ber and angular momentum l a quantum numbers. The 
core states for principal n c and angular momentum l c 
quantum numbers are |c;m) with to £ M. c for M. c = 
{—l c ,...,l c }. Continuum states are approximated as 
free-electron states, i.e., plane waves [see Eq. (JTJ], and 
are denoted by \k) for all k G R 3 flQ8j . Specifically, I 
use momentum-normalized free-electron wave functions 



k) 



1 



i k-r 



(27t) 3 / 2 



(1) 



with the position operator r [108 ] . Taking | k ) for the 
continuum electrons implies that one makes the strong- 
field approximation (SFA) [ill, 0, Ef3], i.e., I neglect 
the impact of the Coulomb potential on continuum elec- 
trons. Only electrons with m € M. for .M = M a fl M c 
are amenable to a core excitation by the x rays in 
electric dipole approximation. The HHG process from 
states | a; to) with m 6 Ma\M, is described within 



the SAE approximation [33j, [34| following Lewenstein et 
al. [ll| because no core electrons couple to these valence 
vacancies in electric dipole approximation. Hence there 
are N/2 = #A4 a + #.M relevant states in the system 
occupied by N electrons. As tunnel ionization of valence 
electrons takes predominant ly place along the linear op- 
tical laser polarization axis [99l llOOj . mostly the valence 
state along the polarization axis is depleted which has 
a magnetic quantum number to = which then couples 
via x rays to core electrons with to = 0. Furthermore, 
the matrix element that is responsible for HHG in the 
Lewenstein model, ( fc | r-e*L | a; to) with the electric dipole 
operator r, vanishes for to ^ and s and p valence states 
from a central potential [112]. There may, however, be 
a weak influence of the valence states with m =/= 0, if 
non-plane-wave continuum s tates are used and electron 
correlations are considered |l!3| . Consequently, I use 
only to = orbitals in the valence and the core in what 
follows and I drop the magnetic quantum number in the 
notation throughout [H, |63| . 

To describe the quantum dynamics of resonant 
x-ray excitation in the course of HHG, one needs to con- 
sider two active electrons: one in the valence or contin- 
uum and one in the core or valence. The spatial one- 
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electron states are occupied by two electrons each with 
opposite spin. To facilitate a core excitation, the spin 
of the valence vacancy and the core electron need to be 
the same as the interaction with the x rays is spin in- 
dependent. Whether the electron with spin up or spin 
down in the valence is ionized does not change the result 
in either case and thus it is sufficient to consider only 
one case [lOJJ. In other words, one needs to consider 
an effective two-electron system formed by one electron 
from the valence and one electron from the core with 
the same spin, i.e., the system has triplet spin. The to- 
tal two-electron state can be factored into a spin part 
and an spatial part. Due to the Pauli exclusion principle 
for indistinguishable electrons and the fact that the spin 
part is symmetric with respect t o elec tron exchange, the 
spatial part is antisymmetric [HI, Il04{ . I use three differ- 
ent classes of spatial two-electron basis states. First, the 
ground state of the two-electron system is given by the 
Slater determinant 

\ac) = -j=[\a)®\c)-\c)®\a)] . (2) 

This is a normalized linear combination of tensorial prod- 
ucts of the two respective one-electron states which is an- 
tisymmetric with respect to electron exchange which in- 
cludes statis tical electron correlations between electrons 
in them [H, Il04j . As core states are energetically well- 
separated from valence states, electron correlations be- 
tween core and valence states are small and neglecting 
them has a minor influence on excitation and ioniza- 
tion e nergies w hich is referred to as core-valence sepa- 
ration jll4l . Ill5j . Second, the valence-ionized states with 
one electron in the continuum for all k £ R 3 and one 
electron in the core state are 

\kc) = ^= [\k)®\c)-\c)®\k)] . (3) 

Third, the core-ionized states with one electron in the 
continuum and one electron in the valence state are 

\ka) = ^= [\k)®\a)-\a)®\k)] . (4) 

The three classes of basis states (0), ©, and dU are 
mutually orthonormal. 

I apply the three assumptions of Lewenstein et al. (llj 
in a somewhat modified way: 

(a) The system is completely described by the three 
classes of states: the ground state |ac), the 
valence-excited states |fcc), and the core-excited 
states | ka). 

(b) The depopulation of the ground state by optical 
laser-induced tunnel ionization and valence ioniza- 
tion by x rays is represented by a decay rate Tq. 

(c) The continuum electron is described as a free par- 
ticle; the influence of the Coulomb potential of 



the residual ion can be omitted. Specifically, also 
the overlap between the free-electron wavefunctions 
and the ground-state state is neglected. 

Further, I need to add the following clauses for the inter- 
action of the x rays with the ion: 

(d) The effect of the x rays — with energies in resonance 
with the level spacing — is confined to the two-level 
system consisting of | k c) and | k a ) in the ion where 
one electron is excited from c — > a or, for ultrain- 
tense x rays, Rabi flops [73, 89- 98] between the two 
levels. The other bound electrons and the contin- 
uum electron are not influenced by the x rays apart 
from (b) and (e). 

(e) Core hole decay and valence ionization by the op- 
tical laser and the x rays are taken into account in 
terms of phenomenological decay rates r a and r c 
for | k c ) and | k a ) , respectively, which are indepen- 
dent of k. 

Point (d) is justified by the fact that the x rays oscil- 
late very rapid ly and cause only a small ponderomotive 
potential [l!2j . Should the second electron be excited 
to the valence vacancy, tunnel ionization is still much 
smaller in (e) than for the neutral atom as the ionization 
potential of the ion is significantly larger than the ion- 
ization potential of the neutral atom. Thus, optical laser 
ionization of the second electron is strongly suppressed 
due to the exponential d ependen ce of the tunnel rate on 
the ionization potential 99, 100]. I treat the interaction 
with the x rays incoherently via decay rates in point (e) 
in contrast to the coherent formulation of the resonant 
x-ray absorption in point (d). 

The phenomenological decay rates of (d) are de- 
termined following Ref. [88l I assume approximately 
monochromatic x rays, i.e., the photoabsorption cross 
sections <7o, a a , and a c for one-photon absorption by 
the atom in the ground, valence-ionized, and cor e-ion ized 
stat e, respectively, can be obtained with Refs. Ill6l and 
Ill7l and are approximately constant over the bandwidth 
of the x-ray pulse for a central angular frequency wx- 
Further, the x-ray intensity is represented by Ix(t) im- 
plying an x-ray photon flux of Jx(t) ~ With the 
optical-lase r-ind uced decay widths Fl.o(£), ^L,a(t) and 
Fh,c(t) [H, llOOt ] . The instantaneous phenomenological 
decay widths of the atom in the three basis states are 
then given by 

r{(t) = r Lii (t) + 9{t) a* J x {t) + S itC 7c , (5) 

for i G {0, a, c} Q with the intrinsic decay width 7 C from 
Aug er and radiative decay of the core-ionized stat e 164 , 
Ill8| and the Heaviside step function 9 and 9(0) = \ [ll9j ]. 
The factor 9(t) emphasizes that the x-ray pulse is zero 
for t < 0. For CW light, I determine decay rates by 
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averaging Eq. ([5]) over an optical laser cycle 



(6) 



The degrees of destruction of the system by the optical 
laser, the x rays, and decay processes up to time t follow 
from Eq. ^ and are quantified by the factor [881 ] : 



Fi(t) = d(t) / n^df 



(7) 



B. Hamiltonian 

The total two-electron Hamiltonian of the atom in two- 
color light (optical laser and x rays) reads 



H — Ha + Hjj + Hx 



(8) 



it consist of three parts: the atomic electronic struc- 
ture Ha, the interaction with the optical laser H^, and 
the interaction with the x rays Hx- 

The one-electron Hamiltonian for the atomic electronic 
structure reads 



h A = \a)e a {a\ + \c)e c {c\ + / \k) — (k \ d 3 fc , (9) 



with the energies e , e c , and ^- of the onc-clcctron states. 
The two-electron atomic electronic structure Hamilto- 
nian follows as 



Ha = h A ' 



r 

1 ® /ja — i | etc) — {ac\ 



(10) 



r 

kc) —{kc\ 



\ka) — {ka\ 



d A k 



with the unity operator for the one-electron Hilbert 
space 1. 

The interaction of one electron with the optical laser 
field is described in electric dipole approximation in 
length form 

L = E h (t) J [|fc)(fc|e L -r»(a|+h.c] d 3 fc (11) 

K 3 

+ E L (t) I f\k) (k\e L -r\k') (k'\ d 3 fcd 3 fc' . 



(12) 



using the electric field of the optical laser 
E L (t) = E 0L cos(w L <) , 



with the peak amplitude _Eol and the optical laser angu- 
lar frequency wl [Z3, HH • The two-electron Hamiltonian 
for the interaction with the optical laser reads 



H 



(13) 



The optical laser couples only the valence electron to the 
continuum because the coupling of the core electron to 
the continuum is negligible. The influence of the optical 
laser on continuum electrons is represented by the last 
term on the right hand side of Eq. (jll|) . 

The interaction of one electron with the x-ray radiation 
is described in electric dipole approximation in length 
form 

h x = Ex(t) [\c)(c\ex-r\a)(a\+h.c] , (14) 
for the x-ray electric field 



Ex(t) = E x cos(w x i) 



Eox 



giwxt + e -iwx*l 



(15) 



with the peak electric field E x and the x-ray angular 
frequency uix [z3i HH ■ The two-electron Hamiltonian for 
the x-ray interaction follows similarly as in Eq. (|13p to 



H x = h x 1 + 1 ® h y 



C. Equations of motion 



(16) 



I make the following ansatz for the two-electron 
wavepacket 



\%t) =a (t)e- iE ° t \ac)+ / [b a (k, t) e~ iE ° 1 \ k c) 



+ b c (k,t)e- i{Ea+ ^ t \ka}] d 3 k 



(17) 



The index on the amplitudes b a (k, t) and 6 c (fc, t) indicates 
which orbital contains the hole and Eq = (ac \ Ha \ac) + 
i ^ = e a + s c . The phase factors in the ansatz (|17[) 
are chosen such that the equations of motion (EOMs) 
are simplified the most, that fast oscillations are sepa- 
rated, and thus the expansion coefficients vary compar- 
atively slowly. I insert \^,t) into the time-dependent 
Schrodinger equation 



ij t \*,t)=H\*,t) 



(18) 



with the Hamiltonian ((8} and project onto the three 
classes of basis states which yields EOMs for the involved 
coefficients. 

First, I obtain the EOM for the ground-state popula- 
tion by projecting onto the ground state ( a c \ yielding 



dt 



oto(t) 



?ao(t)-i / 6„(M) (a\h L \k)d 3 k. (19) 



The second term on the right-hand side can be neglected 
entirely for low ground-state destruction by the optical 
laser [llj . For a higher destruction rate, the second term 
can be neglected as well, if the influence of the optical 
laser is included in Tq [88] . 
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Second, projecting onto (kc\ for all k S R 3 gives the 
EOMs for the recombination of the continuum electron 
with the valence hole which is coupled to core electrons 



I rewrite Eqs. (|20l) and (121]) as a vector equation defin- 
ing the amplitudes 



d_ 

at 



b a (k,t) 



{k 2 -5 + 2I p ~iY a )b a {k 1 t) 



b(k,t) = (b a (k 7 t),b c (k,t)) T 
and the complex symmetric Rabi matrix 



baik'^ik^k') d 3 fc' 



ib c (k,t) (c| /ix.rwa | a) 
io (t) (k\hi,\a) . 



R 



— d — iT a i?ox 
i? x 5 - i r, 



(25) 



(26) 



(20) 



Here, I make the rotating wave approxima- 
tion h x e- iuxt w /ix.rwa Following Ref. HI, 
I introduce the energy ip = — \{s a + e c + lu x ) = —e a + | 
and the detuning of the x-ray photon energy from the 
energy difference between the core and valence lev- 
els S = e a — e c — uix- For zero detuning, Ip corresponds 
to the valence IP of the atom, but is shifted to higher 
(lower) values, if ujx is smaller (larger) than the level 
spacing. 

Third, projecting onto ( k a | for all k £ K. 3 gives the 
EOMs for the recombination of the continuum electron 
with the core hole 



This yields the combined equations 

|L b(k, t) = -~(R+(P + 2I P )l) b(k, t) 
d 



Ehit) W z b{k ' t] 

ia {t)E h (t) p ia 



(27) 



d_ 

dt 



I determine the matrix of eigenvalues of R as 
A = diag(A + , A_) 
with the eigenvalues 

A± = -^(r Q + r c )± M , 

and the complex Rabi frequency (93j : 



b c (k,t) 



-^k 2 + S + 2I p -iT c )b c (k,t) 

i / b c {k',t)(k\h h \k') d 3 k' (21) 



A 4 = 



s + i (r B - r c ) 



7? 2 



(28) 



(29) 



(30) 



ib a (k,t) (a I ftx.RWA |c) 



The matrix of eigenvectors of R is U. With the eigen- 
vectors and eigenvalues, I transform to the basis of x-ray- 
dressed states. I find the new amplitudes as 



In order to evaluate the matrix elements in the EOMs 
that involve continuum electrons, I use the momentum- 
normalized free-electron wave functions of Eq. ([1]). 
Specifically, (k\h L \a) = E L (t)p L ^ a with p L j: a = and the valence ionization fraction 

(fc|e*L • r\a). The optical laser-induced continuum- 
continuum transitions are caused by 



w + 
W- 



= U 



(31) 



(32) 



(k\h lj \k'} = E h (t) (-i)JL6 3 (k'-k), 



(22) 



using the representation 



(23) 



I recast Eq. ((27J) into 

^ /3(fc, t) = ~ (A + (k 2 + 2 I P ) 1) p(k, t) 

+ E h (t) p(k,t) 
-ia {t)E L (t) p L j Q w , 



(33) 



for the Dirac 5 distribution 119]. Further, 



[c\ hx,WNA I a) 



E, 



ox 



R 



■ Px,< 



ox 



(24) 



because U depends neither on k nor on t. 

Equation (|33p is a system of two decoupled partial dif- 
ferential equations. To transform them into two ordi- 
nary differential equations, I introduce the vector poten- 
tial [IT!, Il20| of the optical laser: 



for px.ca = (c|ex • r\a) and the Rabi angular fre- 
quency Rqx [93j. I choose the spatial matrix ele- 
ment px,ca a n d thus i?ox to be real which is feasible 
due to the central potential used to determine atomic 
orbitals fz3l - 



(34) 



Eql . f , _ 
sin(w L t) e L 
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derived from the electric field (IT2|) . The factor e 
ensures convergence of the integral and rj — > + in the 
last equality. As the time dependence of k in Eq. (|33[) 
can be expressed as k! = k — (t ) + Al (f) [ll[ , the sum 
of the left hand side and the derivative with respect to k z 
in Eq. ([33]) represent a total time derivative of /3(fe',t') 
and I have 

A 0(k',t') = - 1 - (A + [(fc - A L (i) + A h (t') f (35) 

+ 2J P ] l) @(k — A\{t) + A^{t'),t') 

- iao(t') £ L (*') pL,jg-A.(t) + ^(t') a ■ 

This is an ordinary first-order differential equation which 
can be integrated exactly [l!9j | yielding 

P±(k,t) = -iw± J a Q (t')E h (t')p hJ ;_ Ih{t)+Ah{t , )a 
o 

-i /(fe-A L (t)+A L (t')) 2 dt" 

xe *' (36) 
xe- i (¥+*0(*- t O dt / j 

by letting t' — t because then k 1 — k and thus f3±(k, t) = 
f3±(k',t'). Next, one can transform back to the ampli- 
tudes of the bare states via b(k, t) = U /3(fc, t). 

D. The time-dependent electric dipole transition 
matrix element 

The emission of light in electric dipole approximation 
is governed by the electric dipole operator which reads 
for one-electron states 

d= J [\a)(a\e n -<r\k)(k\ (37) 

+ \c)(c\e H -f\k)(k\] d 3 fc + h.c. . 

It is defined such that it describes the recombina- 
tion of a continuum electron with the valence and the 
core hole where I omit continuum-continuum transi- 
tions [ll], EH, 13- Here, en is the linear polarization 
vector of the emitted light which is parallel to e^. The 
one-electron electric dipole transition matrix elements 
are PH\afc = H^h • r\k) and p Hc £ = (c|e H • f\k). 
The two-electron electric dipole operator D is defined in 
analogy to Eq. (fT31) by 

D = d(g>i + i(g)d . (38) 

High harmonic emission is determined by the dipole 
transition matrix element between the atomic ground 
state and the valence and core-excited states,^ respec- 
tively. Namely, it is expressed by T>(t) = ( Wo, t\D \ ^ c ,t) 
where D is the two-electron electric dipole operator, 



\^o,t) is the ground-state part of the wavepacket (fTTf 
at time t and \^ c ,t) is the continuum part [13|- Then, 
the time-dependent dipole transition matrix element fol- 
lows from Eqs. (T7]) and (35]) to 

V(t) = (* ,t\D\* c ,t) (39) 
= o5(t) J [p H>a% b a (k,t) + p Hc £ b c (k,t) e- iwxt ] d 3 k 

R3 

= E %(*)• 

?£{+-} 
The dipole components are 

® ij (t) = U ij e- iS ^ u * t a* (t) J Pj^/^M) d 3 k 

K 3 

/ 

= -iC/ t3 u) i e- i ^ Wx 'a*(t) J 'a(i') (40) 

o 

x£ L (t') y pH,p+A L ( t ), P L^ L (*') Q 

R 3 

x e - *^®*'*') d 3 p di' . 

Here, I use the Kronecker delta 5^ c [l!9j | to insert the 
time dependence with wx for the terms involving core- 
excit ed st ates. I also introduced the canonical momen- 
tum jl2l| p = k — A\ (t) and the quasiclassical action 

Sj& t, f ) - \ J(p + Mt")) 2 dt" + + J P ) (t - f ) . 

(41) 

Equation (|40|) can be understood easily term by 
term [TO ]. Here, tUj Pl p+A L (t') a ^' L ^'-' * s ^ e ^ rac ^ on 
of valence ionization at time i' due to tunneling in the 
optical laser field; during propagation of the free electron 
from time t' to t in the optical laser field, the electron ac- 
quires the phase e -1 ^^'*'' the emission of HHG light 
at time t due to recombination of the continuum elec- 
tron with the state i £ {a, c} is governed by p* _ 

the shift of the HHG spectrum by w x for a recombina- 
tion with the core is given by e~ lSi ' c " x *; finally, JTjy is 
the transformation back to the bare-state basis. The ac- 
tion (I4ip differs from the optical laser-only case by 
the summand 4p (t — t') which represents the splitting of 
the core and valence states due to x-ray dressing. 

To calculate T>(t) [Eq. (|39j) ]. I make the saddle point 
approximati on fo r the integration over the canonical mo- 
mentum p [l!9j |. The stationary points of Sj(p,t,t') 
[Eq. (j4Tj) ] are determined from V p Sj(p,t,t') = and 
are given by the difference of the positions of the elec- 
tron at time t and t' [HI . This suggests to introduce the 
new variable t = t — t' , the excursion time of the con- 
tinuum electron. At the stationary points, the canonical 
momentum satisfies 

E 

p st (t,r) = — e L [cos(o; L i) - cos(o; L (t - r))] ; (42) 
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and the quasiclassical action is 

-2^(1-cosKt)) (43) 

-^C(r) cos( WL (2t-r)) , 

with the ponderomotive potential of the optical 
laser U P = jf^ and [l22T| : 

C(r) = sin( WL r) - ^ sin 2 (w L § ) . (44) 

The multiple integral (|4TJj) is then replaced by a single 
integral over the excursion time 

oo 

® ij (t) = -W ij w j e- iS ^"* t a* (t) J V^P^ (45) 

o 

x ao(t -r)Art,i (t,r) e - iSrt ^*' T > dr , 

upon extending the integration over r to infinity and in- 
troducing the atomic terms at the stationary points of 
the quasiclassical action 

A st 4t,r)=E L (t-r)p* HAt{tT)+Mt)i (46) 

X f J L,ft t (t,r)+ J 4 L (t-T)a ' 

The factor \J ^~ 2 ^ '■ stems from the remaining integra- 
tion in the saddle point approximation [l!9j | . The dipole 
components without x rays are obtained from Eq. (|45[) by 
the replacements UijWj e _1<5i cWx * — > 1 and in Eq. (|43l) 
the replacement A± — > —S — iT a . Then the resulting 
equation becomes the same as Eq. (13) of Lewenstein et 
al. apart from a few inessential sign changes, the 
decay width and the additional ground-state amplitudes 
here. 

The atomic terms (|46[) are periodic with the optical 
laser oscillation Tl = %E . anc | can thus be expanded in a 
Fourier series [HI Imf : 

OO 

A stti (t,r)= ]T kiWe- i(2Wi " Kt , (47) 

M= — oo 

with the coefficients 

T L 

b Mli (r) = yI A^Me^+^^'df . (48) 
o 

The bM,i( r ) reflect the properties of the atomic states. 
The definition of the series and the coefficients incor- 
porates the properties of the atomic states: the odd 
and even Fourier coefficients — restricted to an integra- 
tion range from to Tl/2 — can be examined by substi- 
tuting t — » t + 5s Using the properties of Eqs. (IT2"j) . 
(fM]) , (l4*2"j) . and ((11)]) under this replacement I find that 
the argument of the dipole matrix elements is negated. 



To investigate the impact on the atomic dipole matrix 
elements, I express them with Eq. (fTJ and the spatial 
atomic orbital (fi (f) = ( f \ i ) as follows 

Pl,Si = (^372 / e " ife> 3 • r W(0 d 3 r . (49) 

K 3 



The central symmetry of the atomic pote ntial implies 
that the orbitals are parity eigenstates |108j | and a space 
inversion yields p. } — ±p. ^ depending on the parity 
of \ i). Going back to Eq. (|45|) . I find that for i = a, the 
even Fourier coefficients vanish. As a — > c is a dipole 
transition, the states |a) and |c) have opposite parity. 
Therefo re, fo r i = c the odd coefficients are zero. Follow- 
ing Ref. Ill2l . Eq. (|49|) can be decomposed into a sum over 
radial integrals times spherical harmonics times Clebsch- 
Gordan coefficients. 



E. High-order harmonic generation spectrum 

1. Low ground-state depletion 

First, I consider the HHG spectrum for low ground 
state depopulation by tunnel ionization due to the opti- 
cal laser and low valence ionization by the x rays, set- 
ting To = 0. In this case, ao(t) 1 for all times t 
and Eq. (TT9")) can be neglected altogether. As I assume a 
CW optical laser (|12l) , the atomic terms (l4"5"j) are periodic 
in t with period Tl. This implies that the dipole compo- 
nents [Eq. (|45l) ] without the wx-dependent term are also 
periodic with period Tl (in general wx is not a harmonic 
of Consequently, e l5i ' cWx * S)y(t) can be expanded 
into a Fourier series in analogy to Eqs. (|47|) and (|48|). 
The even (odd) Fourier components are present when 
the even (odd) Fourier components of the Fourier series 
of the atomic terms (|4T[1 are nonzero. 

I simplify e i,5i .<=" xt Dij(t) [Eq. (gSJ] by expanding the 
exponential of the quasiclassical action by decomposing 
the i-dependent part of e -lSBt ' J ( t,T ) in terms of the Jacobi- 
Anger expansion 



= J2 (-i) m c- imfl JM, (50) 



with the Bessel functions J m (z) [ll9l |. The general ex- 
pression for the Fourier coefficients of the dipole compo- 
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nents (|45|) without dependence on wx reads: 



in frequency domain by a Fourier transformation 



2) 



ij,2K+8i, a 



- [ e iS ^^ xt ' D lJ (t')e i( - 2K+ ^^^ t ' dt' 



[Jt^e-^) (51) 



f) 

OO 

X f>K--M,i(T) i 

M=— oo 



A/ 



X 



C(r)) e iM «^d7 



where the integration over t' leads to 6k-m,n and I trans- 
lated by M -> K — M. Further, I define 



(52) 



t-2 — ~— (l - cos(cj l r)) 



Inspecting Eqs. (|51|) and (|52l) . I find that there is a 

c 1 2 T dependence of the integrand in Eq. (|5"Tj) . The 

eigenvalues of the Rabi matrix (|29l) have a negative imag- 

i . _.Al t 

inary component — — (r a + r c ) which causes e 1 2 to 

acquire an exponentially decaying factor. In other words, 
the longer the excursion time of a continuum electron, 
the stronger is its impact on the HHG spectrum sup- 
pressed. This effect adds another disadvantage to long 
trajectories versus short trajectories in addition to the 
spreading of the wave packet (|T7|) 0,113 ■ Along the same 
line, multiple returns of the continuum electron to the 
nucleus are diminished. For the case i = a, Eq. (|51l) 
becomes the Lewenstein et al. 11] result [Eq. (18)]— 
apart from constants — for a vanishing x-ray field, because 
then U a j = <5j,+ and Wj — <5j,+. 

The dipole moment T>(t) [Eq. (|39|l ] is obtained from the 
Fourier coefficients of the components (|5"Tj) as follows: 



°(*) = E E [%w e " (2W)utt (53) 
K=o 3 e{+,-} 

+ £> CJ ,2if e- i2J ^ Lt e- iwxt ] . 

The coefficients ([51]) for negative K correspond to the 
inverse process where the HHG photon is emitted before 
the required number of optical laser photons is absorbed 
to satisfy energy conservation [l3j]. The energies of the 
emitted HHG photons are certainly the same as the pho- 
ton energies emitted in the regular process. The inverse 
process is suppressed because the responsible Fourier co- 
efficients have a highly oscillatory integrand and thus are 
very small; they are omitted in Eq. (|55|) . I express T>(t) 



V(t) e lUt dt 



(54) 



2?r E E (®aj,2K+i5((2K+l)u; h -rL) 



k=o je{+,-} 



-D 



cj,2K 



(5(2 K uj l + lo x - SI)) 



usin g the representation (|23[) for the Dirac S distribu- 
tion [HI- 

The HHG spectrum is given by the frequency-resolved 
rate of the HHG photon emission 0, Gj where I con- 
sider here only emission along the x axis with polarization 
vector en = e z , I arrive with Appendix „ at 



dr H (n) 
dQ s 



= 8tt 2 n g (n) w(n) 



with the density of free-photon states 
g(Q) 



(2tt) 3 c 3 



(55) 



(56) 



in the interval [f2; i7+dO] 93] , the solid emission angle fts , 
the speed of light c in vacuum and 



W(Q) = Wi(Q) 



H=0 



Wi(fi) 



E E 

K=0 i£{a,c} 

w"e{+ -} 



(57a) 

\Pi3,2K+5i, a ®if,2K+6 i}a ] 

(57b) 



xS((2K + Si, a ) uj l + 5 itC uj x - fl) 

OO 

£ £ *((2if+l) WL -fi) 

x [®lj,2K+l®cf,2(K-H) +C.C.] . 



(57c) 



The HHG spectrum consist via Eq. (154)) of discrete lines. 
The width of the lines in the HHG spectrum is zero — due 
to the (5 distributions in Eq. (f54|) — despite the widths of 
the singly-ionized states T a and T c . The r a , T c enter 
Eq. (|45p only via t/y, Wj, and Xj which do not depend 
on t and thus cause only a change of the height of the 
lines. The heights are the HH strengths which are deter- 
mined by W(il). Without x rays I obtain 



w(n) 



K=0 



®a+,2K+l | S((2K+l)oj L -n) (58) 



which becomes with Eq. (|55|) — apart from the factor 
of 2 — identical to Eq. (16) in Ref. 0- The experimen- 
tal observable in HHG is the harmonic photon number 
spectrum (HPNS) [13] ■ The HPNS produced by a single 
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atom is obtained from Th(^) by multiplication with the 
duration (square pulse) for which the HPNS is recorded. 
For realistic (Gaussian) optical laser pulses, one needs to 
take into account that the atomic terms can no longer be 
expanded into a Fourier series (|47|) and thus the equa- 
tions in this paper need to be reformulated. 

There are several contributions to the HHG rate in 
Eq. (|55p which I analyze in the following. If wx = 
(2E+l)w L holds for an H G N , then I obtain W(Q) = 
Wi{n) + W 2 , H (ty from Eq. ([57]); the HHG lines due to 
the recombination with the core align with higher orders 
of the recombination with the valence which leads to in- 
terference between both processes mediated by W^ir (£1). 
Otherwise, if wx is not an odd harmonic of ojl, then 
Eq. (dU reduces to W(fl) = Wi(il) and the lines in 
the HHG spectrum belong to two distinct groups: first, 
i — a, HHG light due to recombination with the valence 
vacancy gives rise to lines at the energies (2 K + 1) cjl- 
The x rays induce interference effects in the recombina- 
tion via the two x-ray dressed states with j 6 {+, — }. If 
the x rays are sufficiently intense, the core electron in the 
two- level system may even undergo Rabi flopping (73l . l89l - 
[98j prior recombination with the continuum electron. In 
the bare-states picture, this means that there are mul- 
tiple pathways for the recombination with the valence: 
once, direct recombination and, alternatively, recombi- 
nation after an even number of Rabi cycles of the second 
electron. Second, i = c, HHG light due to recombination 
with the x-ray dressed core vacancy leads to lines at the 
energies 2 K u>l + wx- In other words, the HHG photon 
energies from a recombination with the core are shifted 
by the photon energy ux of the x rays to higher energies. 
In the bare-states picture, an odd number of Rabi cycles 
was performed by the second electron. 



2. X-ray induced ground-state depletion 

Second, I consider the HHG spectrum with ground 
state depopulation due to ionization by the optical laser 
and the x rays by letting To > 0. In doing so, I break 
the periodicity in time imposed by assuming CW optical 
laser and x rays. In other words, the x-ray dressed states 
of Eq. (|33p become an approximation to the x-ray dressed 
states of the nonperiodic system. In the ensuing Sec. IIII1 
I will examine the HHG process for arbitrary x-ray pulses 
and overcome this limitation. This will facilitate to assess 
the results of this subsubsection and study the impact of 
a finite x-ray pulse. 

To determine the time-dependent ground-state popu- 
lation, I need to solve Eq. (|19l) where I neglect the sec- 
ond term on the right-hand side and include the destruc- 
tion by tunnel- ionization in Tq. For a constant ionization 
rate To starting at t = 0, the solution for the ground-state 
amplitude reads 



a (t) = 6(-t)+< 



Est 

2 



0(t). 



(59) 



x-ray dressed states can only be used for t > 0. Before 
the x-ray pulse, the optical laser causes the usual HHG 
for valence electrons. The duration of the x-ray pulse is 
Tp. I will focus on the part with the x rays and disre- 
gard the additional contribution from t < and t > Tp 
to the valence-recombination part of the HHG spectra. 
Due to the lack of temporal periodicity, the HHG spec- 
trum cannot be expanded into a Fourier series any longer 
and I need to ca lculate the Fourier transform of Eq. (|45l) 
instead [H [H|: 

r P 

%(Q) = [® l3 (t)e int dt 



11,,'r, j v /np!e- i ^M (60) 
o 

N=-oo 
oo 

M-— oo 

where I expanded the exponential of the quasiclassical 
action at the stationary point and translated by M — ¥ 
M — N. I define the line shapes 



-i [(2 (M+N)+S iia ) w L +5 iiC u x -!l] t 



x a* (t) a(t - t) At . (61) 

Equation (|60p has a straightforward interpretation. 
Comparing it with the result without ground-state deple- 
tion (|54[) . I see that the 8 distribution line shapes in the 
latter equation are replaced with finite-width line shapes 
in the former. Ground-state depletion thus "adds" a fi- 
nite width to the Fourier coefficients. 

Using the solution (f59j) for the ground-state depletion, 
I find for the line shape ([ST]) : 



(0,r) 



e^(l 



-Tq— i (^M,N,i— Tp 



r + i(n 



M.NA 



fi) 



(62) 



Clearly, a(t) is not periodic in time. Additionally, the 



Neglecting the factor e~ T , I see that /ijw,iv,i(^) T ) peaks 
at Q = (2 (M + N) + 5i.a) wl + <5i )C wx. In other words, 
the peaks are at the positions of the harmonics with- 
out ground-state depletion; for i = c, the harmonics are 
shifted by with respect to the harmonics for i — a. 
Est T 

The factor e 2 leads to a deformation and shifting of 
the peaks. Above all, hM,N,i(S^-: T ) h as a finite width; 
its real part is proportional to a Lorentzian of a width 
of4r . 

The HPNS for a single atom — the probability to find 
a photon with specified energy — along the x axis follows 
from the Appendix and is given by 

dP H (^) 



An s An 



(63) 
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with X>(fi) = [compare with Eq. (|55|) ]. 

i£{a,c} 

3. Cutoff law 

The cutoff of the HHG spectrum with x-ray dressing 
is changed quite significantly with respect to the opti- 
cal laser-only case. First, inspecting Eq. (|39|) . I see that 
the terms for the recombination of the continuum elec- 
tron with a core hole are shifted by u>x towards higher 
energy compared with the terms for the recombination 
with a valence hole. Neglecting the changes in the atomic 
terms (|4"6")l between the two cases, I see that also the cut- 
off of the HHG spectrum resulting from the core hole 
increases by u>x- 

Second, the phenomenological semiclassical cutoff law 
of an optical laser-only HHG spectrum [H, [ll[ is given by 

n cut = 3.17C/p+X P , (64) 

with the valence ionization potential Zp . X-ray dressing 
introduces, however, a changed situation. Namely, I ob- 
serve that Eq. contains two terms each for valence 
and core recombination with j — ±. Associated with 
the x-ray dressed states j • = ± are the ionization poten- 
tials -#+ip. In other words, from Eq. (|64|) I see that the 
two terms of valence or core recombination have different 
cutoffs; the higher of the two determines the cutoff of the 
HHG spectrum. Consequently, there can be a few har- 
monics close to the cutoff that stem only from the term 
with the higher ionization potential. 

III. HIGH-ORDER HARMONIC GENERATION 
WITH ARBITRARY X-RAY PULSES 

I would like to examine HHG in the presence of a gen- 
eral x-ray pulse as they are obtained from x-ray free elec- 
tron lasers which are frequently operate by the SASE 
principle [68l470l | using the PCM [7ll473l |. The formalism 
of the previous Sec. HT1 relies on CW x rays [Eq. (Tbfj) ] by 
going into a dressed-states picture [Eq. (l33|) ] which is no 
longer practicable for arbitrary x-ray pulses. In this sec- 
tion, I pursue a derivation that circumvents the dressed- 
states picture. Expressions which have counterparts in 
Sec. Illl have the same symbol with a prime attached. 

A. Hamiltonian, wave function ansatz, and 
equations of motion 

In order to describe the interaction with arbitrary 
x-ray pulses, I need to replace the interaction with 
CW x rays h x [Eq. 03}] by 

h' x = 6(t)(e x -r)^^- [ e iKi+ w W) +e -i(,x^xW)] _ 

(65) 



Here, E' ox (t) is the time-dependent amplitude of the 
pulse start ing a t t = and fx(t) is the time-dependent 
phase [7l|, 1 123( | . The two-electron Hamiltonian for the 
x-ray interaction H x has the same form as Eq. (|16p 
with hx replaced by h' x . 

The newly arisen dependence on the time-dependent 
phase of the x-ray field (|65l) needs to be taken into ac- 
count in the ansatz for the two-electron wavepacket (|17l) 
which is modified by multiplying the phase factor of | k a ) 
additionally by e _i Vx(t) . With this ansatz I obtain EOMs 
which are very similar to the ones obtained previously 
apart from the decay widths T' (t), T' a (t), and T' c (t) which 
are now time dependent. First, the ground-state popula- 
tion is formally identical to Eq. (|19p . Second, also the 
recombination of the continuum electron with the va- 
lence hole is formally identical to Eq. (|20p . Third, in 
the recombination of the continuum electron with the 
core hole, Eq. ([2"T]) . the first term on the right hand 
side needs to be augmented by the summand — 2tfix(t) 

with ip x (t) = ^r 1 - The modified Eqs. (J20|) and (ED 
are recast in terms of a matrix equation formally identical 
to Eq. (g7]) yielding 

f t b'{k,t) = - X - (R'(t) + (k 2 + 2/p) 1) b'(k,t) (66) 

+ Mt) m- b'{k, t) - io'(t) E L (t) Ph>ia Q . 

with b'(k, t) = (,?,;-' I) where the instantaneous Rabi ma- 

v ' \b' c (k,t)' 

trix [compare with Eq. (|26p ] reads 

R'(t) = -idiag(r:(i),r^) - 2i¥*(t)) ( 67 ) 
- 5er z + R' ox (t)cr x , 

with the Pauli matrices <r x and er z [92|, Il08| 

and the instantaneous Rabi frequency i?Q X (i) = 
2 (c|^ x ,rwaI«> = Kx(t) Px, ca for h' x e- ("* '+**(*)) w 
rwa = (^x • r) — ^— using the rotating wave approxi- 
mation [93(. 

B. Integration of the equations of motion 

In order to solve the system of EOMs (|66p , I need to 
realized that, in contrast to the EOMs for CW x rays 
[Eqs. (H9l) and ([27]l]. the Rabi matrix R'(t) [Eq. ^} is 
time dependent. A decoupling of the two recombination 
amplitudes, as in Eq. (|31l) is, therefore, no longer possible 
by going to an x-ray dressed-states picture because the 
eigenvector matrix of R'(t) does not commute with the 
time derivative. However, I can decouple the EOMs (|66p 
in a somewhat different way. I begin by removing the 
diagonal elements of R (t) by the substitution 

b' a {k,t)=b'i{k,t)^ st - L ^ L (68) 
b' c (k,t) = %(k,t)e-5 st - 2 , 
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with b"(k,t) = Inserting the ansatz (|68|) into 

Eq. (JBBJ leads to 

5 f) = (-A i^ x (f) E x (i) - i (y + /p) l) ?'(*, t) 

(69) 



0t 



9 



+ E L (t)—b"(k,t) 



ia'(t)E L (t) p u%a 



i ... Fq(t) 

e -2 <>*+— J— 



The matrix 



e A W 
-A(t) o 



(70) 



with the phase-amplitude difference 

A(i) = -i<5 t + i + ^ (M*) - r c (*)) (71) 

in Eq. (|69p does not depend on momentum but does de- 
pend on time. Therefore, a diagonalization of S a (t) leads 
to time-independent eigenvalues 

A' = diag(A' + ,A'_) =tr z , (72) 

and time-dependent eigenvectors 

/ e A(t) _ e A(t)> 



1 



1 



(73) 



I transform to the eigenbasis of E^i) — in analogy to the 
dressed-state basis used before [Eq. ([33)) ] — yielding 

1 0'(k, t) = (if -\t) U'(t) - ~ R' ox (t) A' 



-i^ + Jpjlj^fc,*) (74) 

+ E L (t)-^-J'(k,t) 

-ia'(t)E h (t) p Ljfc ->'W ■ 

I added U'- X (t) U{t)f3'{k,t) on both sides of Eq. fffg]) 
where the new amplitudes [compare with Eq. pip ] are 



0'(k,t) = (^j^j) = U'-\t)b"{k,t) . (75) 

I use a dot over symbols to denote a time derivative ' = 
-4: ; the transformation rest reads 

u'- 1 (t)U'(t) = ^(-l + cr x ), (76) 
using Eq. (|7ip and the valence ionization fraction 
w'(t) 



tfM _ (w'+{t) 



w'_(t) 



u'-\t) 



Fajt) 







(77) 



= I e2 -5t-ivx(t) + £i 2 ^ ^ 1 



[compare with Eq. (|32l) ]. 

In order to decouple the two differential equations (|74l) , 
I need to approximate the transformation rest (f76|) by 



its diagonal elements tf l {t) U'(t) ps — 1. I refer 
to this as zeroth order and attach a superscript W to 
the involved coeffic ients . The resulting equations can be 
integrated exactly [119] analogous to Eqs. (|35l) and (|36l 
yielding 



/4 (0) (M) = i / a'W{t')E h {t')p u% _ A , t 



Here, the quasiclassical action [compare with Eq. (|41|)] is 



5j(fc,t,t') = i / (fc-I L (t) + i' L (i")) 2 dt" 



+ y (e(t)-e(t')) 

-A(A(t)-A(f)) 

+ / P (t - o , 

and pulse area [11] of the x rays is 

/ 

e(t) = 0(t) / ^x^Odf 



(79) 



(80) 



C. Time-dependent dipole moment and high-order 
harmonic generation spectrum 

The time-dependent two-electron dipole transition 
matrix element follows from Eq. (I38[) by transform- 
ing (3'(k,t) [Eq. (|T51) ] back to the bare-state ampli- 
tudes b'(k,t) with the inverse of Eqs. (|fi"51) and (I75"j) : 

2y(i) = (* / 0J *l^l*c.*> = E ( 81 ) 

2^{a,c} 
J6{+,-} 

where I denote by I^q,*) the ground-state part of the 
modified wavepacket (fTT)) and by \^' c ,t) the continuum 
part [lj]. 



i. Zeroth order 

I introduce the dipole components 2)^ (i) using the 
zero-order coefficients 

^° ) (t) = Ul j (t)e^) a '(0)* {t) J Pu . n ^) { k,t)d s k, 

(82) 
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and the time-dependent phase 

<f>i(t) = (-1)*" 1 1 - Si, c ux t + i . (83) 

The dipole components (|82l) are not periodic in time 

and thus the (t) cannot be expanded into a Fourier 
series. Inste ad, I need to take the Fourier transform 
of them Il23l |. I introduce the canonical momen- 
tum p — k — A^(t) and the excursion time r = t — t' , 
expand the atomic terms into a Fourier series [Eqs. (|46|) . 
(|47p , and (|48p] and make the saddle point approxima- 
tion [119] which yields 

OO i 

= - /7fc*2!.-»sw 



OO 

x E i JV ^(^C(r))e iJV '^ (84) 

N=-oo 

OO 

x X] ^/-jv^t) h' M iJ (il, t) dr , 

M=-oo 

where I expanded the exponential of the quasiclassical 
action at the stationary point, used Fq(t) from Eq. ([52]) . 
C(r) from Eq. (|gg]l, translated by M -4- M - N, and 
defined 



x e -|Ai(e(t)-e(t-r)) e i0i(t) ^ 85 ^ 

xe -i(A W -AM) fl ,(0), (() 

x a'^-^l^^-T) df . 

As in Sec. HIE 21 I disregard the contributions to the 
HHG spectrum from valence recombination in Eq. (|85[) 
before the x-ray pulse sets in at t = 0, i.e., I assume a op- 
tical laser pulse which starts at zero time. As I assume 
a finite x-ray pulse, I restrict the integration (I85[) to a 
optical laser pulse duration of Tp (square pulse). This 
result (|84p has a close resemblance to the dipole compo- 
nents for CW x rays with ground-state depletion 
However, the eigenvectors U (t) transform from a differ- 
ent basis to bare states. 

To find an expression for Eq. (pjj) , I need to solve the 
modified Eq. (p~9|) where I neglect the second term on the 
right-hand side because I assume a low tunnel ioniza- 
tion rate. For a x-ray induced ionization by an arbitrary 
pulse starting at t = 0, the solution for the ground-state 
amplitude reads 

o'(°)(t) = 0(-t)+e- 2 6{t) , (86) 



(3J, dZ3, I®, and (US]): 



, ,, <icij _ 7" Fo(t)+Fo(t-T) 
h'M,N,i,j(n,T) = ( ' 1} 2 - y e 3 



x e -i [(2 (M+N)+S l . a )ui L +S i , c wx-O] « 
x e 5 (px(t)-<?x(t— r)) e -i<5«,c fx(t) (g7) 
-FaW-F.W+r^i- r)+r c (t— r) 

x e 4 

x e -^i (ew-e(t-r)) ^ 



Ignoring the impact of the other factors in the equation, 
the line shape peaks at il = (2 (M+N)+5i <a ) wl+4cWx- 
The other factors differ substantially from the result (|62l) . 
The reason for these contributions to occur is the only 
partial diagonalization of the instantaneous Rabi ma- 
trix (157)) . The suppression of the lineshapes is caused 
by the destruction of the ground state at the time of 
tunnel ionization Fo(t — r) and at the time of recombi- 
nation F o(t) of the continuum electron with the ground- 
state part. Further, the term depending on F a (t), 
F c (t) represents the destruction of the intermediate hole 
state by x-ray absorption or Auger decay in between tun- 
nel ionization and recombination. Inspecting Eq. (|29[) . 
I see that via Eq. (JTJ) this term and the term depend- 
ing on the pulse area (l80|) can be understood to result 
t 

from J dt where I retain only the leading order in 

t-r 

the expansion of the complex Rabi frequency X± with 
respect to i?ox which is taken to be large with respect 
to all other parameters in /it, i.e., /j, ss i?ox- Finally, the 
probability of the HHG photon emission along the x axis 
follows in zeroth order from Eq. (I63[) . 



2. First order 



I have obtained the solution of the coupled equa- 
tions (|66p taking only the diagonal elements of the trans- 
formation rest (|76l) into account. Higher order approxi- 
mations to the amplitudes /?'(fc, t) can be obtained by it- 
eration, i.e., the zero-order result obtained in Sec. MI (JTI 
is used in Eq. ([74"]) to specify the offdiagonal elements of 
Eq. ((75)) . For the first-order correction f3'^(k,t) to the 
wave function in the expansion 



similarly to Eq. (|59l) . This expression is inserted into 
Eq. ((551) to obtain the lineshape. I find with Eqs. (I72"j) . 



(? (jfe, t) « /3'(°) (jfe,t) + (*,*), (88) 
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I find from Eq. ([74 



|^(M) = (-i(^ + / P -iA(t))i 



--R' 0X (t)A') (3'(k,t) 
+ E L (t)^0'(k,t) (89) 

-ia'W(t)E L {t)p LiSa vf(t) 
+ |A(i) a x ^ \k,t) . 

These equations are still decoupled and can be solved 
along the lines presented in this section [Eq. ([75]) ] by 
inserting the expansion (|88[) which yields 



= / |A(t / )e- 1 ^±(* At 'Jo- B ^°)(jfe,t / )dt / • 
o 

(90) 

Next, the zero-order solution [Eq. (|78|) ] is inserted into 
Eq. (HO), the two actions are united and treated with 
the saddle point method [l!9j . Defining two excursion 
times t' = t — t' and r = t — t" leads with the expan- 
sion of the atomic terms (|47l) to the first order correction 
of the dipole components (|82p . (|84p that manifests in a 
first-order correction of the line widths (|85p in Eq. (154") ; 
the expression of the first-order correction of the dipole 
components reads 

©^(fi) = -i / UUt)J*®a'®>(t)tr* A ® 



xe -2^ e (*) ^ K M ,s,iM (91) 

JW,iV=-oo 
x e -i [(2 (M+JV)+5i,a) WL-H] t ^ 



with the integrals 



Mt-r') pi A,e(t-Q 



(92) 



and 



xi%(^C(r)) e iAr ^ (93) 
x a'^(i-r) to' -(t-r) 



x e2 



A(t-r) e -i Xj6(t-r) dr 



The first-order correction (I9ip differs significantly from 
the zero-order expression (|85[) . Differences arise due 
to the additional integration over r'. The dependence 



of A(t — t') and 0(t — r') on r' accounts for the phase- 
amplitude change due to the offdiagonal matrix elements 
in Eq. ([76]) . The form of the integral JM,N,i,j(t) has a 
certain resemblance to the integration over t in Eq. (|84p ; 
it represents the effective amplitude accumulated for all 
tunneling instances at t — r in the allowed interval [0;t]. 
The integral K,M,N,i,j(t) quantifies the impact of the 
change of the phase and amplitude via the offdiagonal 
elements onto the dipole components. 

I obtain closed-form expressions for 2)^^(0) by using 
Eqs. (72} , (73), (77), (SD), (S3), and (SB) which yields 

t p 

J.— 

"ij w ' 2 



© 1) (r») = -i^A- / e- 2 4 



x e2 



x e 



o 

e~2 



M,N=-oa 
-i [(2 (M+A r )+5 i-a ) w L +5i. c wx-n] < 



(94) 



X )Cu,N,i,j(t) dt , 



with the integrals 



1 



^M,N,i,j(t) = J 2 [-i5 + i^x(t-r') 
o 



+ f«('-^')-fc(t-r') ] e iA,e(t-r') (Q5) 

x \jM,Ns,j{t) — Jm,n,ij(t')] dr' , 



and 





xi%(|C(r)) e iw ^ r 



(96) 



Fo(t-T) Fa ( t -T) + Fc(t-T) _i_ 

e 2 + 4 e 2 

e -|A 3 .e(t-r) d 



tpx(t— t) 



The probability of the HHG photon emission along the 
x axis follows by inserting the sum of the zero- and first- 
order corrections into Eq. (63). 



IV. CONCLUSION 

In this paper, I discuss theoretically the impact of 
x-ray excitation on the HHG spectrum. I derive equa- 
tions for a two-electron model that is based on the one- 
electron model for HHG by Lewenstein et al. [llj. The 
first electron from the atomic valence tunnels through 
the potential barrier into the continuum where it prop- 
agates freely. Additionally, I consider a second electron 
from the atomic core which is driven by intense x rays 
that are tuned to the core-valence resonance in the ionic 
remnant which is formed by tunnel ionization. For high 
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x-ray intensities, this electron Rabi flops between the va- 
lence hole and the core. The returning first electron thus 
sees a superposition of valence-hole and core-hole states 
and recombines with the valence or the core hole emit- 
ting HHG radiation that is characteristic of this super- 
position. I analyze the HHG spectrum and identify the 
parts due to recombination with the valence and the core. 
The latter light is shifted by the x-ray photon energy to 
higher energies. 

In my treatment, I focused on the HHG spectrum of 
a single atom. However, in experiments, a macroscopic 
sample is used to study HHG. The propagation of the 
optical laser, xuv, and x-ray radiation through the gas 
transforms the spectrum significantly due to interferences 
effects caused by coherent emission of light by the sam- 
ple. The phase matching [22j shall be examined in future 
studies specifically under the objective of the impact of 
SASE x rays on the HHG process; estimates of the yield 
of x-ray-boosted HHG show that a considerable output 
can be achieved from a macroscopic medium [88j. 

The model of Lewenstein et al. [ll| describes excel- 
lently qualitatively HHG spectra but tends to overesti- 
mate the H HG y ield significantly by roughly two orders of 
magnitude [l24j | . These shortcomings of the model with 
respect to quantitative accuracy have been suggested to 
be m itiga ted by a few methods, using the Ehrcnfcst theo- 
rem [124{ for the reco mbin ation step, the eikonal-Volkov 
approximation (EVA) [1251 ] , and considering electron cor- 
relations [13, Ell- In the case of a SAE, a numerical inte- 
gration of the time-dependent Schrodinger equation [l26j | 
is a route to obtain an accurate HH yield. However, 
in our case, a numerical integration of a two-electron 
Schrodinger equation woul d be required with substan- 
tially increased complexity |l26l |. 

My prediction offers novel prospects for nonlinear 
x-ray physics which is practicable unlike simultaneous 
multi-x-r ay- p hoton absorption that exhibits tiny cross 
sections jl27l |. Namely, the HHG process determines a 
time window in which x-ray interaction may occur. Rabi 
flopping is a highly nonlinear fundamental process and 
becomes experimentally feasible with x rays for the first 
time using FELs 73, 94-98]. The occurrence of Rabi os- 
cillations may have a much clearer signature if combined 
with the HHG scheme compared with the case if only 
x rays are considered. HHG and x rays may also repre- 
sent an optical gatin g to perform frequency resolved op- 
tical gating (FROG) [l28| with x-ray pulses thus offering 
the long-sought after puls e characterization for chaotic 
SASE FEL light [HI |l30| but requires further theoreti- 
cal research. The nove l sche me also makes single attosec- 
ond x-ray pulses [l3ll 1 1321 ] and attosecond x-ray pulse 
trains [35[ feasible using the same methods that are used 
with conventional HHG (88[. Above all tomographic 
imaging of core orbitals with HHG comes into reach |25l - 
I271 . 129I . l30j . Finally, HHG has been used to generate fre- 
quency combs in the XUV; potentially, such HHG spectra 
can be boosted by x rays in order to extend frequency- 
comb-based spectroscopy to the x-ray regime [l33l 1 134| 



provided that one has identically-shaped x-ray pulses 
which are carefully synchronized to the optical laser. 
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Appendix: Harmonic photon number spectrum 

1. Quantum-electrodynamic formulation 

The recombination of the continuum electron with the 
hole in the cation leads to the emission of a HHG photon. 
This fluorescence cannot be treated with semiclassical 
fields; instead, quan tum electrodynamics is required flU 
EE El Ell HI EHJ. The quantized electric field of the 
HHG light is in dipole approximation 

4 = (-i) £ ^ 4 - & . (A- 1 ) 

with the photon energy D,^ of the qth mode, the normal- 
ization volume V used to quantize the electromagnetic 
field, the polarization vector e*H along the z axis, and 
annihilation a - and creation at operators for photons in 

mode q. The Hamiltonian for the HHG fluorescence Hn 
is constructed based on Eq. (|3"5)l by replacing the polar- 
ization vector with En. The free HH field is 

where I set the energy of the vacuum to zero. The quan- 
tum electrodynamic Hamiltonian for the description of 
fluorescence reads 

Hqed = Ha + Hem + #h , (A-3) 

with the atomic Hamiltonian Ha from Eq. (|10p . 

The coefficients a(t), b a (k,t), and b c (k,t) from the 
EOMs of Sec. Ill CI were obtained with no HHG photons 
present and are assumed to be known here. To determine 
the amplitude of fluorescence, I need a similar wavepacket 
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to Eq. (|17l) which includes the photon number states of 
the HH held: 

=^c f (t)e- i ^ t |ac)®|l,-) + |* c ,t)®|0) . 

(A.4) 

Here, | $ Cl t) is the continuum part of Eq. (|17[) . the vac- 
uum state is 1 ) and 1 1^) = at j ) . 

2. No ground-state depletion 

This new wavepacket is inserted into the time- 
dependent Schrodinger equation and projected on the 
state (ac \ <E) (l<f| where no ground-state depletion is as- 
sumed, i.e., To = 0. This leads to an EOM for the am- 
plitude to find a fluorescence photon similar to Eq. ([TO]) : 

-± c&) = n ? c q ~(t) - i V{t) , (A.5) 

with the the time-dependent two-electron dipole tran- 
sition m atrix element (|39|) . The equation is integrated 
directly jl!9| yielding 

t 

ctft) = -y/^Sk e - in ^ J e in **' £>(*') dt' . (A.6) 
-t 

The expression for T>(t') can be replaced by a Fourier 
series expansion (|53|) . In the limit t — > oo, the time 
integral in Eq. (|A.6j) becomes a Fourier transform f|54[). 
I carry out the time integration for finite t using 

t 

J e~ ixt ' dt' = 2tsmc(xt) , (A.7) 
-1 

with the sinus cardinalis for z G C: 

(l ;z = 

sincz = < sinz . (A. 8) 

I— 5 ^°- 

The amplitudes are 

. oo 

j£{+-}K=0 

x [&aj,2K+i sine ((2 K + l)uj L t- Q,$t) 

The probability to measure a photon is given 

by \ cg(t) \ 2 . For large t, I transform it with the rela- 
tion [jog : 

lim t sinc 2 (xi) = tt<j"(x) . (A.10) 

t— >oo 

The integral over x of t sine (x t) sine ((x + y) i) for y ^ 
vanishes in the limit £ — ^ oo. Therefore, the cross terms in 



the probability vanish and are dropped in the following. 
The rate of photon emission is given by d (2t) l c <rWI 2 Us- 
ing the fact that free-photon states are spaced densely, 
I arrive at the rate in Eq. (|55p upon integrating over 
V times the density of free-photon states g(Q) [Eq. (|56l) ] 
and collecting the results in the energy dependent rate 
in Eq. (I55|) . The extra factor of 2 in Eq. (|55|) accounts 
for the two possible spin states in the valence orbital 
with to = 0. The HPNS per atom is obtained from the 
rate by multiplication with the optical laser pulse dura- 
tion for a square pulse [see discussion after Eq. (|55|) ]. 



3. Ground-state depletion by x-ray absorption 

The situation of the previous subsection changes as 
soon as I consider ground-state depletion by setting To > 
0. In this case, the HHG emission is not periodic any- 
more. Instead of a rate, I, therefore, calculate the HPNS 
for a single atom directly, i.e., the probability for an atom 
to emit a HHG photon. 

To describe the probability amplitude of HHG fluo- 
rescence photons, I make a product ansatz for the co- 
efficient c^t) = ^Ijy in the wavepacket in Eq. (|A.4[) 
where a(t) is the ground-state amplitude (|5§1) [or Eq. (|8"6"j) 
for a'(t)]. Replacing T a by -T in H A from Eq. (JTOJ) to 
account for loss of norm in T>(t), the division by a(t) coun- 
ters the increase of norm due to —To such that C<f(i) is not 
affected. This leads to an EOM in analogy to Eq. (|A.5[) 
which has the solution 

a|_.^/.^*«'" (a.ii) 

-t 

x a*{t')- 1 V{t')dt l . 

The time integration over the decay width in the expo- 
nent yields two factors; the first factor simply is a(i) _1 
and cancels with the factor on the left-hand side. The 
second factor forms the dipole matrix element with the 
remainder of the integral such that I arrive at 

c q <t) « -yf^ e-^ 4 • (A.12) 

The integral in Eq. (|A.11[) becomes the Fourier trans- 
form of V(t') in the limit t -> oo . The HPNS from a 
single atom follows from Eq. (|A.12[) by taking the modu- 
lus squared, letting t — > oo, and replacing the mode q by 
V g(fl) [Eq. dSSD] yielding Eq. @3) [the extra factor of 2 
stems from the two spin states]. In Sec. V of Ref. [TJ an 
average HHG-emission rate is constructed from | "Z7 ( O ^*) | 2 
for nonvanishing ground-state depletion by multiplying 
its value at the positions of the harmonics with the area 
of the harmonic peaks. 
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